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Abstract 

We study responses of the Brans-Dicke field due to gravitational collapses of scalar field pulses using 
numerical simulations. Double-null formalism is employed to implement the numerical simulations. 

If we supply a scalar field pulse, it will asymptotically form a black hole via dynamical interactions 
of the Brans-Dicke field. Hence, we can observe the responses of the Brans-Dicke field by two different 
regions. First, we observe the late time behaviors after the gravitational collapse, which include 
formations of a singularity and an apparent horizon. Second, we observe the fully dynamical behaviors 
during the gravitational collapse and view the energy-momentum tensor components. 

For the late time behaviors, if the Brans-Dicke coupling is greater (or smaller) than —1.5, the 
Brans-Dicke field decreases (or increases) during the gravitational collapse. Since the Brans-Dicke 
field should be relaxed to the asymptotic value with the elapse of time, the final apparent horizon 
becomes time- like (or space- like). 

For the dynamical behaviors, we observed the energy-momentum tensors around ui ~ —1.5. If 
the Brans-Dicke coupling is greater than —1.5, the T uu component can be negative at the outside of 
the black hole. This can allow an instantaneous inflating region during the gravitational collapse. If 
the Brans-Dicke coupling is less than —1.5, the oscillation of the T vv component allows the apparent 
horizon to shrink. This allows a combination that violates weak cosmic censorship. 

Finally, we discuss the implications of the violation of the null energy condition and weak cosmic 
censorship. 
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1 Introduction 



Einstein gravity is the most successful theory to describe gravitation. However, if we include quantum 
effects, Einstein gravity cannot be a fundamental theory, and hence a quantum theory of gravity is 
required. This quantum gravity should be approximated by the Einstein gravity in the low energy limit, 
but before approaching the Einstein gravity, there may be some effects that modify the Einstein gravity. 
In fact, a theory of gravity does not definitively have to be Einstein gravity. There can be many variations 
and we choose one theory via observations. We then need to find possible modified gravity theories and 
check their possible implications. If quantum gravity is approximated by a modified gravity theory, and 
if the theory has experimental implications, it will be profoundly important to understand the nature of 
quantum gravity. 

One of the most canonical modified gravity theory candidates is the Brans-Dicke theory [T] : 



where R is the Ricci scalar, <f> is the Brans-Dicke field, and cu is the Brans-Dicke coupling constant. His- 
torically, the Brans-Dicke theory was suggested to implement Mach's principle. The roll of the gravitation 
constant G is replaced by the Brans-Dicke field 1/0. Since the Brans-Dicke field can be dynamic, the 
strength of gravity can be relatively different for each points. It makes that the equivalence principle 
has more freedom than the original Einstein gravity, since the equivalence principle holds for each points 
independently. Then, intuitively, we can recover the Einstein gravity in the large |w| limit, since the field 
dynamics of the Brans-Dicke field will be highly suppressed in the limit. 

Moreover, the Brans-Dicke gravity has other theoretical motivations. Notably, it appears that it is 
the simplest form of modified gravity using a scalar field. Therefore, some important intuitions of the 
Brans-Dicke theory can be applied to other possible modified gravity theories. For example, we may apply 
similar intuitions to dilaton gravity, non-minimally coupled theory, etc. 

The Brans-Dicke theory can be motivated by string theory inspired models. For example, the Brans- 
Dicke theory is the weak coupling limit of dilaton gravity (lu = —1), where the dilaton field is a direct 
consequence of string theory [2]. Also, the Brans-Dicke theory is a weak field limit of the Randall-Sundrum 
model [3J, where, on the positive tension brane, lu is sufficiently large; on the negative tension brane, 
lu > —1.5 [4 5 . Of course, the observed limitation of the coupling is u > 4 x 10 4 6 ; however, it is still 
meaningful to study other values of to since they can be realized in the context of string theory, and hence 
they can show properties which are allowed by string theory. For further discussions, see Appendix A. 

If the Brans-Dicke field is added to Einstein gravity, it can work as a source of an exotic matter. In 
the context of cosmology, there have been some discussions that the Brans-Dicke field can be a candidate 
for dark matter or dark energy [7]. In particular, it can be useful to study exotic matters that violates the 




(1) 
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null energy condition. The theory can allow some strange geometry, e.g., wormholes [5]. Therefore, if we 
know how to deal with the Brans-Dicke field so as to obtain a sufficient amount of exotic matters, then 
we can reasonably make an assumption of exotic matter, since the exotic matter may thus be naturally 
obtained from a well-known and not artificial field. 

Although the study of the Brans-Dicke field is important, solutions to the Brans-Dicke theory are more 
difficult than those involving Einstein gravity. Of course, we can find some analytic solutions. However, 
to observe the real dynamics when a black hole is formed and evolved, more detailed calculations are 
needed. The collapsing matter will affect the background, and the background will affect the Brans-Dicke 
field. The affected Brans-Dicke field will then give a back-reaction to the background geometry. These 
processes are not easy to understand in an analytic way and it is necessary to use numerical calculations 
for their elucidation. 

Previous papers have studied Brans-Dicke black holes using ADM formalism 9 . However, in this 
paper, we use double-null formalism [l0][Tl][12] for the numerical calculations [l3][T4], thus making it 
easier to study the causal structures and also to plot all possible fields including the energy-momentum 
tensors. 

The purpose of this paper is to study responses of the Brans-Dicke field during the gravitational 
collapse of a matter. We prepare a flat background with a pulse of a normal scalar matter field and we 
then observe the response of the Brans-Dicke field. We thereupon observe that dynamical responses of 
the Brans-Dicke field are highly non-trivial. 

This paper is organized as follows: in Section [2J we discuss the details of our numerical simulations; in 
Section [31 we discuss the results of our simulations and report some interesting observations; in Section 21 
we document observations in terms of two issues (late time behaviors and dynamical behaviors) and discuss 
their physical interpretations; and in Section [5j we summarize the possible causal structures and discuss 
the physical meanings and implications. 

2 Model for Brans-Dicke theory 
2.1 Brans-Dicke theory 

The Lagrangian of the Brans-Dicke theory with a scalar field becomes [T] 



where R is the Ricci scalar, is the Brans-Dicke field, and $ is a minimally coupled scalar field with 
potential V($). Here, u is the Brans-Dicke coupling constant which is a free parameter of the theory. 




(2) 
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The Einstein equation becomes as follows: 



G\xv — 87rT^, + 8tt— ^- = 87rT M „, (3) 



where the Brans-Dicke part of the energy-momentum tensors are 



and the matter part of the energy-momentum tensors are 



The field equations are as follows: 



%^S""-V'(i) = 0, (7) 



where 



T = T p. (8) 

2.2 Implementation in double-null formalism 

We use the double-null coordinates 

ds 2 = -a 2 (u, v)dudv + r 2 (u, v)dVL 2 , (9) 

assuming spherical symmetry, where u is the retarded time, v is the advanced time, and 9 and tp are 
angular coordinates. 

We follow the notation of [H] [10] [11] [12] : the metric function a, the radial function r, the Brans-Dicke 
field (f>, and a scalar field S = \^4tt<&, and define 

ct a 

h=—, d=—, / = r u , g = r v , w = 6 U , z = <j>,v, W = S, u , Z = S v . (10) 
a a 

The Einstein tensors are then given as follows: 

Gun = --(f, u -2fh), (11) 

r 

Guv = ^ (4r/,„ + a 2 + 4/ 5 ) , (12) 

G vv = --(g, v -2gd), (13) 
r 

Gee = + v) • (14) 
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Also, we can obtain the energy-momentum tensors for the Brans-Dicke field part and the scalar field 
part: 



-iBD 



rpBD 



rpBT) 



rriBD 

1 ee 



1 



z,u _ gw + fz 



w 



-Z.u + 



2 / ko> 2 4> ' 4ira 2 (j) 



(gw + fz) - 



uj r 
a 2 



wz. 



(15) 
(16) 
(17) 
(18) 



UU 



UV 



VV 



-—w z 

4tt0 ' 



4tt0 ' 



(19) 
(20) 
(21) 
(22) 



4> 2ira 2 (j) 

To implement double-null formalism into the numerical scheme, it is convenient to represent all equa- 
tions as first order differential equations. Note that 
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The Einstein equations for a. uv , r„, and the field equation for <j> are then coupled: 
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where 



21 = 
23 = 



1 



£ = — /z — gw 



r 2 (j) 00 2r 4> 
a 2 fg ( 4-7rr 
4r r 

4?rr /WZ 



20 2 



+ — TZ-^(gw + fz), 



3 + 2w V 2tt 
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(25) 
(26) 
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After solving these coupled equations, we can write all equations: 



t , , r „ ru 



/,„ - 2fh-—(w, u -2hw)--W 2 -—w 2 , (28) 



g, v = 2gd-—{z^-2dz)--Z 2 -—z 2 , (29) 
fg a 2 gw + fz u 2n (WZ 2 \ WZ 

d > u = h '* = ^ + ^ + — W wz + TzTWA~^~ )~~ ' m 

fg a 2 gw + fz 2nr (WZ 2 \ litre? 

!U = J ' = -T-^-^2^ + j3T^{-^- 2aV ) + — v > (31) 

If „ 47rr /WZ „ 2t 

/ i £ _ . ■ / ^^ ,-^1 



= -^ 5W + /Z+ 3T^^- 2 ^))' (32) 



including the scalar field equation 



f Z 9W 



Z U = W V = -- ira^V'(S). (33) 

r r 

Now, equations of a. uv , r jU „, <p tUV , and S. uv parts can be represented by first order differential equations. 
We can then implement the same integration scheme used in previous papers [11] [12] to solve the Brans- 
Dicke theory. We used the second order Runge-Kutta method [16] . Tests of the convergence are provided 
in Appendix B. 



2.3 Initial conditions and free parameters 

We need initial conditions for all functions (a, h, d, r, f, g, <f>, w, z, S, W, Z) on the initial u — u\ and v = V{ 
surfaces, where we set u\ = Vi = 0. 

We have gauge freedom to choose the initial r function. Although all constant u and v lines are null, 
there remains freedom to choose the distances between these null lines. Here, we choose r(0,0) = ro, 
/(w, 0) = r u o, and g(0, v) = r„o, where r u o < and r„o > such that the radial function for an in-going 
observer decreases and that for an out-going observer increases. 

First, we assume that the gravitation constant G = 1/0 is asymptotically 1. Then, <fi(u, 0) = 4>(0, v) = 1 
and w(u, 0) = z(0, v) = 0. 

We use a shell-shaped scalar field, and hence its interior is not affected by the shell. Thus, we can 
simply choose S(u, 0) = A and a(u, 0) = 1. Also, W(u, 0) = h(u, 0) = holds. Then, since the asymptotic 
mass function |17] 

m(u,v)^ r -{l + 4 r -^), (34) 

should vanish at u — v — 0, it is convenient to choose r„o = — 1/2 and r v o = 1/2. 

We need more information to determine d,g,z, and Z on the v — surface. We obtain d from 
Equation ([3H|) . g from Equation ([3"Tj) . z from Equation (|3"2"|) . and Z from Equation ([3"3"]h respectively. 
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We can choose an arbitrary function for S(0, v) to induce a collapsing pulse. In this paper, we use 

(35) 





A . . 




K) 













for < v < Vf and otherwise S(0, v) — 0, where Vf is the width of the pulse and A is the amplitude. Then 
we obtain 

Z(0,v) = -— sin 2 (n-) (36) 

Vf \ Vf J 

for < v < Vf and otherwise Z(0,v) = 0. This implements one pulse of energy (T vv ~ Z 2 ) along the 
out-going null direction by a differentiable function Z{Q,v). 

Also, from Equation (|29l) . we can use d = rZ 2 /2g(f> on the u — surface and we obtain d(0, v). By 
integrating d along v, we have a(0,?;). 

We need more information for h, /, w, and on the u = surface. We obtain /i from Equation (|3T)]) . 
/ from Equation (|3"Tj) . w from Equation ([3^1l . and W from Equation (133)) . respectively. This finishes the 
assignments of initial conditions. 

We choose ro = 10, leaving the three free parameters (u),A,Vf), where u> is the Brans-Dicke coupling 
parameter, A is the amplitude of a pulse of the scalar field, and Vf is the width of the pulse. Here, we 
assume that there is no potential term in the matter field side. 



3 Responses of the Brans-Dicke field due to gravitational col- 
lapses 

3.1 Simulation results 

We report the simulation results in Figure [TJ We fixed A = 0.75 and Vf — 20, and varied oj = 
10, 1, 0, — 1, —1.4, —1.6, —2, —10. Note that if lj = —1.5, because of Equation ©, the equation becomes 
singular. 

Figure [1] contains contour diagrams of the radial function r for each parameters. We plotted singular- 
ities, trapping (apparent) horizons, and Cauchy horizons. 

• If the radial function becomes 0, since the r = curve is space-like and all equations become singular, 
it is reasonable to interpret this region as a singularity. 

• To define a black hole using local geometry, we may use apparent horizons [15] [15] or trapping 
horizons [20j for an out-going observer, i.e., r „ = 0. If there is inflation at some point in space-time, 
the inflating region can be defined by an in-going observer, i.e., r , u = 0, so that the in-going observer 
sees an increase of the area function. We plotted two horizons: r „ = and r u = 0. 
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Figure 1: Contour diagrams of the radial function r. A = 0.75, Vf — 20, and u = 
10, 1, 0, —1, —1.4, —1.6, —2, —10. Spacing of each contour is 1. Here, we plotted singularities, horizons 
(r >v = and r u — 0), and Cauchy horizons. 
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Figure 2: The integration domain of the double-null simulation. Region (A) is for the late time behavior; 
(B) is for the dynamical behavior. 

• If an apparent horizon contracts to a singularity, we can find a null cutoff line where we cannot 
calculate, because we do not have sufficient information to determine this region via the singularity; 
hence, we can identify the region as a Cauchy horizon |18]|19j. 

3.2 Observations 

In this subsection, we list some interesting observations for Figure [1] In the next section, we interpret the 
observations. 

1. Eventually a black hole is formed, i.e., a space-like singularity and an apparent horizon are observed. 

2. If | a; | is sufficiently large, the causal structures become similar to each other and the apparent 
horizons are always space-like. See ui = 10 and —10 cases. 

3. As uj decreases, if uj > —1.5, the apparent horizon approaches the event horizon in the time- like 
direction. See u = 1, 0, —1, —1.4. 

4. If uj < —1.5 and u) is sufficiently small, the apparent horizon oscillates and eventually approaches 
the event horizon in the space-like direction. See u) = —1.6, —2. 

5. If w > —1.5, the causal structure can contain an instantaneously inflating region, i.e., there can exist 
an r_ u = horizon. As uj approaches —1.5, the r u = horizon appears from upper u to lower u. 
Compare uj = — 1 and uj = —1.6. 
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6. If lj < —1.5, the apparent horizon can shrink to a singularity and form a Cauchy horizon. However, 
eventually a black hole will be formed. See uj = —1.6. 

Observations 1, 2, 3, and 4 pertain to asymptotic behaviors after the collapse of matter. We call them 
late time behaviors. Observations 5 and 6 pertain to fully dynamical behaviors during the gravitational 
collapse. We call them dynamical behaviors (Figure [5]). We discuss and interpret the late time behaviors 
and dynamical behaviors in the following section. 

4 Interpretations 

4.1 Late time behaviors and dynamics of the Brans-Dicke field 

From the work of previous authors using Brans-Dicke theory, we know that a gravitational collapse should 
lead to the emergence of an Einstein black hole asymptotically [3T] . As the black holes emerge to Einstein 
black holes, some interesting observations can be made. To understand these observations, in this section, 
we discuss the large |w| case and the uj ~ —1.5 case independently. 

4.1.1 The large \uj\ limit 

If |w| is sufficiently large, the field equation of the Brans-Dicke field (Equation ([5])) becomes a free scalar 
field equation. Then, if <f> was asymptotically 1 and initially had no dynamics, it will not be affected by 
gravitation, since the field equation is a free scalar field equation. (See the cases of uj = 10 and uj = — 10 
in Figure [3]) Therefore, the dynamics during the collapse of matter will be similar to that of Einstein 
gravity and we can ignore the effect of the Brans-Dicke field. This explains why the black holes in the 
large |w| limit resemble Einstein black holes. 

Note that, in these cases, space-like apparent horizons are observed. If there is no supply of matter, 
the horizon should be null. However, during the matter collapse, there is some scattered energy along 
the out-going direction, and part of this energy will come to the black hole later. This explains why the 
horizons are still space-like and not null after Vf — 20. 

4.1.2 uj near —1.5 

If uj is near —1.5, we will see dynamical effects of the Brans-Dicke field. Figure [3] shows that if uj is greater 
than —1.5, it decreases during the matter collapse; if uj is less than —1.5, it increases during the matter 
collapse. We discuss and interpret these phenomena below. 
We have the equation for the Brans-Dicke field: 




(37) 
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A=0.7a, vi=20 






Figure 3: Plot of the Brans-Dicke field </>. A = 0.75, v { = 20, and w = 1 0, f , 0, -f , -f .4, -f .6, -2, -fO. 
If the color changes, it shows that the Brans-Dicke field has dynamical behaviors. If uj > —1.5, then 
decreases during the matter collapse; if uj < —1.5, then <fi increases during the matter collapse. 
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Figure 4: Plot of the scalar field S. A — 0.75, Vf = 20, and to = — 1 and w = —1.6. These figures show 
that, during collapse of matter, S decreases for both the in-going and out-going directions. Therefore, 
S.uS. v > holds during matter supply. For other parameters, the behaviors are similar during the matter 
collapse. 

During the matter collapse, the S^ u S yV term is dominant. Also, in our initial condition, S^ U S^ V is always 
positive during the matter collapse (Figure [4]). Then, the total sign of <p^ uv is determined by the sign of 
3 + 2oj. Therefore, if u> > —1.5, the total sign of cf>, uv becomes negative, and hence the Brans-Dickc field 
tends to decrease. If ui < —1.5, by the same reasoning, the Brans-Dicke field tends to increase. 

After the matter supply ends, the Brans-Dicke field will be relaxed to the asymptotic value, since the 
black hole should approach an Einstein black hole. Hence, if ui > —1.5, the decreased Brans-Dicke field 
tends to increase for an out-going observer around the horizon. Note that the area function in the Einstein 
frame is <j>A [12], where A is the area of a black hole in the Jordan frame. We know that, according to 
the area theorem, the area of the Einstein frame (<pA) will always increase and be space-like. Then, 

5{<j>A) > 0. (38) 

If there is no collapse of additional matter, 8(<f>A) ~ and <f>A ~ const. Therefore, if increases along 
the horizon, the area should shrink and be time-like. If u) < —1.5, by the same reasoning, we can explain 
that the horizon contracts along a space-like direction. 

4.2 Dynamical behaviors and energy conditions 

Now we turn to a discussion of fully dynamical behaviors during the matter collapse. We have observed 
interesting behaviors of r t = or r u = horizons. These behaviors should be consistent with energy- 
momentum tensors. Therefore, in this section, we look at details of the energy-momentum tensors and 
energy conditions that cause the dynamical behaviors. We discuss two different cases, u> > —1.5 and 
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u> < -1.5. 

To understand the dynamical behaviors, we plot the energy-momentum tensor components (Figure [5]). 
For any observer whose four velocity is n M , the observer feels the energy density by 

T^n". (39) 

Therefore, if an observer moves along an in-going null direction, the observer feels the energy density by 
T uu ; if an observer moves along an out- going null direction, the observer feels the energy density by T vv . 
Note that the null energy condition is violated if T uu or T vv has a negative region. 
Around horizons (r u = or r, v = 0), we can simplify Equations (|28[) and (|29[) by 

r, uu = -4irrT uu , (40) 
r, vv = -AirrT vv . (41) 

We know that asymptotically r „ < and r „ > 0; i.e., initially, an in-going observer sees a decrease of 
area and an out-going observer sees an increase of area. Using these observations, we remark on two cases 
where the null energy condition should be violated. 

• We can define the outer part of an r u = horizon if the sign of r M changes from — to + along an 
in-going null direction; the inner part of an r u = changes its sign from + to — along an in-going 
null direction. Then, around the outer r u = horizon, r uu should be positive, and hence T uu < 0. 

• If r vv < around the (outer) r v = horizon, the horizon is space-like; if r jVV > around the r ; „ = 
horizon, the horizon is time-like. Therefore, to see a time-like r„ = horizon, T vv < is required. 

Note that, T uv component is related to the vacuum energy or cosmological constant term. Therefore, 
if T uv > 0, then it may cause an increase of the area; if T uv < 0, then it may cause a decrease of the area. 
We can write Equation (|3"TT) around the r u = horizon as 

r tUV =-^-+4<7rrT uv . (42) 
4r 

If there is inflation, an out-going null observer sees an increase of area. This requires that r jU „ > 0, and 
hence T uv should be sufficiently large: 

< T uv . (43) 



167rr 2 



If the inflation tends to end, the opposite situation should occur and T uv should be sufficiently small so 
that r yUV < 0. Therefore, we need to check the sign of T uv components if there is an r „ = horizon. 
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Figure 5: The energy-momentum tensors for (lu — —1.4, A — 0.75, Vf = 20). 



«=-1.4, A=0.75, vr=20 
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Figure 6: The energy-momentum tensors of the matter sectors (T® u /(j>, T® v /4>) and the Brans-Dicke sectors 
ftuu i T vv D ) for (w - -1.4, A = 0.75, «t = 20). 
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4.2.1 w>-1.5 

In FigurejSJ we can see the energy-momentum tensor components of the u> > —1.5 case. The outer r. u — 
horizon was possible since the T uu < condition holds. Also, r „ = could bend in the time-like direction 
since the T vv < condition holds. Around the r u — horizon, the T uv component changes its sign. We 
can interpret that, as the area increases, some effective vacuum energy can be gained from the Brans-Dickc 
field. However, as the effective vacuum energy changes its sign, the area ceases to increase and eventually 
decreases for an in-going observer. We call this an instantaneous inflation driven by gravitational collapse. 
We can compare the energy-momentum tensors by two sectors: 

T» v = ?f+T™, (44) 

where T*^ and T^ D are defined by Equations [5] and H] In Figure [IjJ we compare these two sectors. During 
the matter supply, the matter part and the Brans-Dicke part are equally important for T vv . However, in 
other cases, the Brans-Dicke part dominates. 

4.2.2 w<-1.5 

In Figure [3 we can see the energy-momentum tensor components of the u> < —1.5 case. One interesting 
observation is that T vv has a negative region when the r^ v = horizon shrinks. The negative part of 
T vv originates from the negative part of (Figure [5]). As v increases, the sign of is changed by 
+, — , +. This is an opposite behavior than the u> > —1.5 case; in this case, was changed as — , +, — . 

These differences are consistent with the late time behaviors. When the matter collapse finishes, T^j 3 
dominates the total energy-momentum tensor. Hence, in the to < —1.5 case, the apparent horizon is 
space-like; in the us > —1.5 case, the apparent horizon is time-like. 

4.2.3 Violation of weak cosmic censorship 

Finally, we discuss the violation of weak cosmic censorship [T5] for the ui < —1.5 case. The existence of a 
Cauchy horizon that is not hidden by an apparent horizon implies the violation of weak cosmic censorship. 
However, it is interesting to check whether the existence of the Cauchy horizon can be maintained infinitely 
or there should be a singularity and an apparent horizon, as in the (w = —1.6, A = 0.75, Vf = 20) case. We 
argue that if the asymptotic location of the singularity has an upper u coordinate relative to the location 
of the Cauchy horizon, wc cannot see the singularity forever. 

We briefly remark on a property of our double-null coordinates. 

• Eventually, our black holes should approach a Schwarzschild black hole, and in this limit, we can 
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use the Kruskal coordinates [33] so that 

uv = (r-2M)e X p(J^), (45) 

where M is the mass of the Schwarzschild black hole. Then, the singularity is uv — — 2Af and the 
horizon is uv = 0. 

If we choose u = at the horizon, then for a fixed v, the distance between the horizon and the 
singularity in terms of the u coordinate becomes Au oc l/v. Then, as v increases, the singularity 
will approach the horizon in terms of the u coordinate. 

• The location of the apparent horizon r „ = will follow the case of a Schwarzschild black hole. Also, 
we know that the horizon of a Schwarzschild black hole is only determined by the mass. Therefore, 
if we fix A and Vf , and hence we fix the asymptotic mass, for any u, we will see asymptotically the 
same u coordinate for the apparent horizon. 

• Therefore, the final location of a singularity is the same if we fix A and Vf. This is observed in 
Figure [T] 

Therefore, first, to find the location of the singularity, for given A and V{, we will see a black hole of 
the u> = — 1.4 case. Q Second, we compare the location of a Cauchy horizon for the to = — 1.6 case with the 
same A and same V{ . If the location of a Cauchy horizon in lj = —1.6 has smaller u than the location of a 
singularity in lj = — 1.4, this confirms that the existence of a Cauchy horizon can be maintained infinitely. 

In Figures [51 and [TU1 we present contour diagrams with various initial conditions. 

In Figure IH1 we fixed V[ — 20 and changed A = 0.55,0.65,0.85,0.95. Also, we compared uj = —1.4 
and w = —1.6. In the cases of A = 0.85 and 0.95, the locations of the singularity between u = —1.4 
and u> = —1.6 converge to the same location. This confirms our previous remark on the double-null 
coordinates. Therefore, in the cases of A = 0.55 and 0.65, we can check that the singularity will not 
appear below the Cauchy horizon, and hence the Cauchy horizon will be maintained infinitely. 

In Figure [TUJ we fixed A = 0.75 and changed Vf = 15,30,40. Also, we compared lu — —1.4 and 
lj = —1.6. As in the previous paragraph, we can see that V{ = 30 and 40 cases allow infinite extension of 
the Cauchy horizon. 
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Figure 9: Contour diagrams of the radial function r. «f = 20, A = 0.55, 0.65, 0.85, 0.95, and u> = —1.4 and 
— 1.6. Spacing is 1. 
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Figure 10: Contour diagrams of the radial function r. A = 0.75, Vf = 15, 30, 40, and u = —1.4 and 
Spacing is 1. 
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Figure 11: Classification of solutions. 
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5 Conclusion 



5.1 Summary of results 

We study responses of the Brans-Dicke field due to gravitational collapses of scalar field pulses using 
numerical simulations. In this paper, we used the double-null coordinates to implement the numerical 
simulations. 

If the absolute value of the Brans-Dicke coupling constant |w| is sufficiently large, the responses of the 
Brans-Dicke field are negligible and eventually an Einstein black hole will be formed. However, if uj is 
sufficiently small such that the Brans-Dicke field becomes sufficiently dynamic, we can see fully dynamical 
back-reactions of the Brans-Dicke field. 

If we supply a scalar field pulse, it will asymptotically form a black hole via dynamical interactions of 
the Brans-Dicke field. Hence, we observed the responses of the Brans-Dicke field by two different regions 
(Figure [5]). First, we observed the late time behaviors after the gravitational collapse, which includes 
formations of a singularity and an apparent horizon. Second, we observed the fully dynamical behaviors 
during the gravitational collapse and viewed the energy-momentum tensor components. 

For the late time behaviors, if uj is greater (or smaller) than —1.5, the Brans-Dicke field decreases (or 
increases) during the gravitational collapse (Figure \S§ . Since the Brans-Dicke field should be relaxed to 
the asymptotic value with the elapse of time, the final apparent horizon becomes time-like (or space-like) 
(Figured]). 

For the dynamical behaviors, we observed the energy-momentum tensors around uj ~ —1.5. If uj > 
— 1.5, the T uu component can be negative at the outside of the black hole. This allows an instantaneous 
inflating region during the gravitational collapse (Figure [3]). If uj < —1.5, the oscillation of the T vv 
component allows the apparent horizon to shrink (Figure [7]). In this case, we could find a combination of 
initial conditions that allows a violation of weak cosmic censorship (Figures and [TU]) . 

In Figure [TT1 we summarize and classify all solutions in this paper. 

(A) If \uj\ is sufficiently large, the r^ v — horizon is space-like. 

(B) If uj > —1.5 is sufficiently small, then the r v = horizon can be time-like after the matter collapse. 

(C), (D) If uj > —1.5, there can be r u = horizons. After the matter collapse, the r >v = horizons are 
time-like. 

1 Oi course, we could compare other ui values. However, the lo = —1.4 case clearly shows asymptotic time-like horizons. 
The location of a singularity should have upper u coordinate relative to that of a time-like horizon. Therefore, the ui = — 1.4 
case is useful to clarify the location of the singularity. 
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(E) If uj < —1.5, we can find an initial condition that the r„ = horizon shrinks to the singularity and 
eventually a space-like r„ = horizon and a singularity appears. We can see a Cauchy horizon. 

(F) If uj < —1.5, we also find an initial condition that the r, v = horizon shrinks but not to a singularity. 
Then, we see an oscillation of the r v = horizon. 

(G) In uj < —1.5 case, we can find an initial condition that the Cauchy horizon can be maintained 
infinitely. 

5.2 Discussion 

Finally, we remark on some physical implications of our results. 

First, in the uj > —1.5 limit, we could see a violation of the null energy condition for the T uu component. 
This condition is a necessary condition to induce a bubble universe using a false vacuum bubble |12) . It is 
necessary to check further calculations to ascertain whether this is in fact useful; but one lesson is that, in 
Brans-Dicke theory or similar type of theory, we may apply the back-reaction of a gravitational collapse 
to induce the negative T uu part. 

Second, in the uj < —1.5 limit, we could see a violation of weak cosmic censorship. Even though the 
Cauchy horizon extends infinitely, it is reasonable to speculate that a Schwarzschild black hole will be 
formed beyond the Cauchy horizon. However, it is also true that we cannot determine what happens 
beyond the Cauchy horizon. These situations were not observed using numerical simulations by previous 
authors [9]- Therefore, there should be further studies for this clear counterexample of weak cosmic 
censorship. 
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Appendix A. Physical realizations for various u 

For observational tests, it is known that the value of uj should be greater than 4 x 10 4 [6]. However, in 
various physical situations, small uj parameters can be allowed. Even though the small uj is not for our 
universe, if small uj is allowed in fundamental theory and may be realizable in somewhere of multiverse, 
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and if such small value of lu has some implications, e.g., violation of unitarity or cosmic censorship, the 
study of responses of the Brans-Dickc field for small lu will have theoretical importance. 

One example is the dilaton gravity which has the effective action as the following form [2] : 

S = 1^=1 I d d+1 x^—ge-* (R + (V$) 2 ) , (46) 



2A 

where d is the space dimensions, A s is the length scale of string units, R is the Ricci scalar, and $ is the 
dilaton field. If we define <\> as 

(47) 



Xt 1 %nG d+1 

where Gd+i is the d + 1 dimensional gravitation constant, then we obtain the Brans-Dicke theory with 
lu = — 1 limit. If there are higher loop corrections from string theory, there will be other coupling terms of 
<f> and hence the Brans-Dicke theory should be modified. Hence, the correspondence between the dilaton 
gravity and the Brans-Dicke theory is only for weak coupling limit (e - * 3> 1). However, it is reasonable 
to think that Brans-Dicke theory with lu = — 1 limit is a good toy model to study phenomena of dilaton 
gravity. 

In the first model of Randall and Sundrum [3] , they introduced two branes to explain the hierarchy 
problem. Because of the warp factor between two branes, we obtain a positive tension brane and a negative 
tension brane in an anti de Sitter background. According to Garriga and Tanaka [4], each branes can be 
described by Brans-Dicke theory in the weak field limit with the lu parameter 

3 



ui 



| (e ±s/l - l) , (48) 



where s is the location of the negative tension brane along the fifth dimension, I = y— 6/A is the length 
scale of the anti de Sitter space, and the sign ± denotes the sign of the tension. To explain the hierarchy 
problem, we require s/l ~ 35. Then we obtain sufficiently large lu on the positive tension brane and 
lu > —3/2 on the negative tension brane. However, in principle, s/l can be chosen arbitrarily and hence 
one may guess that various lu near —3/2 may be allowed by models from brane world. 
For more discussions, see [5]. 



Appendix B. Convergence and consistency tests 

In this appendix, we report on convergence and consistency tests for our simulations. We used lu = — 1, 
A = 0.75, and V{ = 20 case. 

For convergence, we compared finer simulations: 1x1,2x2, and 4x4 times finer. In Figure IT2l we 
see that the difference between the 1 x I and 2x2 times finer cases is 4 times the difference between the 
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V 



Figure 12: Plots of errors with different step sizes. Here, we plot |r(i) — ?"(2) IA(2) and 4|r( 2 ) — r '(4) l/ r *(4) 
along a few constant u lines, where r( n ) is calculated in an n x n times finer simulation than rny Since 
two curves are folded for each u slice, we confirm that our simulation converges to second order. 



2x2 and 4x4 times finer cases, and thus our simulation converges to second order. For u = 10 and 
u = 15, the curves appear to diverge, since there is a singularity. However, even before the singularity, 
the error is < 0.1%. 

We also compared constraint equations (Equations (|2"51) and (|2"9")l ) using 



2 , ru w 2 



\U + \2fh\ + |MK„| + \2hw\) + \$w*\ + \^w*\ °' (49) 



2<j,\\\" J ,u\ 1 i 10" I i | 2 02 

|g,, - 2gd + ^(z,„ - 2ciz) + r.Z 2 + j^z 2 
| fli0 | + |2 5 d| + |M|z,„| + \2dz\) + \%Z*\ + 



0. (50) 

Figures [T3] and [TJ] show the constraint equations. Around the singularity, the constraint equations increase, 
but still the value is less than ~ 1%. For Figure [TJ] the constraint equation initially appears to be large, 
because initially both the numerator and denominator, in Equation (|50[) . are too small (less than ~ 10~ 10 ). 
Therefore, these plots show that our simulations hold the constraint equations sufficiently. 
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